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Abstract 

We review the phase structure of a three-dimensional, non-compact AbeHan 
gauge theory (QED3) as a function of the number N of 4-component massless 
fermions. There is a critical Nc up to which there is dynamical fermion mass 
generation and an associated global symmetry breaking. We discuss various 
approaches to the determination of Nc, which lead to estimates ranging from 
Nc — I to Nc = 4. This theory with N — 2 has been employed as an effective 
continuum theory for the 2D quantum antiferromagnet where the observed Neel 
ordering corresponds to dynamical fermion mass generation. Thus the value 
of Nc is of some physical interest. We also consider the phase structure of the 
model with a finite gauge boson mass (the Abelian Higgs model). 
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1 Introduction 



We review the dynamical generation of fermion mass in a three-dimensional, non- 
compact Abelian gauge theory (QED3) and its generalization to the case of a finite 
gauge boson mass. Attention was first drawn to this problem by Pisarski [1]. Some 
years ago, a study concluded that with the gauge symmetry unbroken, mass genera- 
tion will occur, corresponding to the breaking of a certain global "chiral" symmetry, 
if and only if the number N of four-component fermions is no larger than a criti- 
cal value Nc, estimated to be 4 [2, 3]. Since then, most studies have agreed that 
there is indeed such an Nc, but there has been little agreement as to its magnitude. 
Estimates have ranged from that of Refs. 2 and 3 to as low as A?c = 1 [4]. This 
uncertainty is typical of our rudimentary understanding of most strongly coupled 
quantum field theories. 

The value of Nc is a potentially important real-world question. It has been 
suggested that QED3 with N = 2 massless fermions can provide a continuum de- 
scription of a 2D quantum antiferromagnet [5, 6]. Marston [7], Laughlin and Zou 
[8], and others have noted that the observed Neel ordering corresponds to dynamical 
fermion mass generation and its associated chiral symmetry breaking in QED3 with 
N = 2. More recent discussions of these ideas are provided by Kim and Lee [9], and 
Kleinert, Nogueira, and Sudbo [10]. For another application of non-compact QED3 
to condensed matter systems see Ref. [11]. The suggestion that (non-compact) 
QED3 provides a continuum description of the 2D quantum antiferromagnet has, 
however, been challenged in a recent series of papers by Herbut et.al [12]. 

Several approaches have been brought to bear on the determination of Nc 
in QED3. The original study [2, 3] was based on a questionable use of the 1/N 
expansion along with a continuum gap equation. Continuing to non-integral values 
of N, it led to the estimate Nc ~ 128/37r^ ~ 4.3. This approach also led to some 
speculation on the nature of the phase transition as a function of N [13], namely 
that because of the long-range force it is of infinite order. A lower value of Nc 
is suggested by a conjectured inequality based on the counting of thermodynamic 
degrees of freedom [4]. It indicates that Nc can be no greater than 3/2. It is likely 
that the determination of Nc will be settled only by numerical lattice studies. These 
studies have so far suggested that Nc could be as low as 2 [14, 15], but a conclusive 
determination will need more refined simulations. 

We first review the features of QED3 and discuss its symmetry properties and 
behavior for large N where the 1/A^ expansion may be reliably employed. We discuss 
fermion mass generation, which takes place at low values of N. The conjectured 
constraint of Ref. 4 is then described and applied to QED3. We review recent 
lattice computations, discuss their sensitivity to finite size effects, and to provide 
some perspective on these simulations include a rough numerical estimate using the 
continuum gap equation. We also discuss dynamical fermion mass generation in a 
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related theory, the Abehan Higgs model in three-dimensions with the same fermion 
content as QED3, but where the gauge field is massive. We then summarize, revisit 
the relevance of these models to condensed matter systems, and describe some open 
questions. 



2 QED3 and its Symmetry Breaking 

The Lagrangian of the model is 

^ - 1 
L = Y.^j{ip)^j-jF^uF'''^, (1) 

where = id^j, + eA^ is the covariant derivative, e is the gauge coupling, and V'j is 
a set of N 4-component fermion fields. To explore the phase structure as a function 
of N , it is convenient to keep fixed the dimcnsionful quantity a = Being supcr- 
renormalizable, this theory is UV-complete, rapidly damped at momentum scales 
beyond a. As an effective theory at lower momenta, one could nevertheless take the 
view that higher momentum physics is unknown, being integrated out into a tower 
of higher-dimension operators associated with the scale a. We do not include such 
terms , which could modify the critical value of N for chiral symmetry breaking. We 
return to this question in the summary when we discuss possible condensed matter 
applications of QED3. 

For large AT, the theory remains weakly coupled at all momentum scales. 
This can be seen by computing the gauge boson propagator in the large-A limit 
and extracting from it an effective, dimensionless running coupling a{k). The large- 
N form of the Euclidean propagator is {k"^ -|-a/c)^^, where k is the magnitude of the 
Euclidean three-momentum and the second term arises from the N fermion loops. 
Multiplying the propagator by the coupling and by one power of k to make a 
dimensionless quantity, we have 

_ _ e^k 8 a 

This expression exhibits asymptotic freedom at large momentum, and conformal 
symmetry with an IR fixed point of strength 8/A as A; ^ 0. For large N, the 
coupling is always weak, and no dynamical fermion mass generation is expected. 

For finite N, on the other hand, the infrared coupling becomes strong, gov- 
erned by a strong IR fixed point, and fermion mass generation becomes possible. 
In the absence of fermion mass, the global symmetry associated with the N four- 
component fermions is U{2N). The dynamical generation of an equal mass for all 
the fermion flavors would break the global symmetry to U (N) U(N) giving rise 
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to 2A^2 Goldstone bosons [1]. More generally, mass terms for the Dirac spinors can 
violate parity (P) and time reversal symmetry (T), and the gauge field admits a cor- 
responding P- and T-violating Chern-Simons mass term. These are not generated 
spontaneously [16]. We concentrate here on the parity conserving case. 

Several studies over the years have led to the conclusion that a parity con- 
serving mass is indeed generated dynamically for this theory providing that A'^ is 
no larger than a critical value Nc [2]. An upper limit is expected, since, as we have 
noted, the theory is weakly coupled at all momentum scales in the large-N limit. 
Since the Lagrangian of Eq. 1 is UV-complete, damped rapidly at momentum scales 
above a, the value of Nc is completely determined by the conformal, IR behavior. 

3 A Conjectured Constraint 

A conjectured constraint on the infrared structure of asymptotically free gauge 
theories [4] can be utilized to analyze the phase structure of QED3. This constraint, 
which takes the form of an inequality, states that for a wide class of such theories, 
the number of IR " degrees of freedom" , defined using the thermal free energy per 
unit volume F{T) of the theory, is less than or equal to the number of corresponding 
UV degrees of freedom. 

These degree-of-freedom counts are defined by the quantities fm and fuv, 
given in terms of F{T) by taking the the zero- and infinite-T limits respectively 
of the quantity — ^^^ fid). Here f{d) is a function of the number of space time 
dimensions d, defined such that the contribution from a free bosonic field is 1. The 
inequality is then 

fiR < fuv. (3) 

In general, fm and fuv will exist providing that the theory is governed by IR and 
UV fixed points. The conjectured inequality was restricted, however, to asymptoti- 
cally free theories, in which fuv can be computed from free field theory. It has been 
applied to a variety of four-dimensional theories, and found to be satisfied whenever 
fjR can be reliably computed [4]. This includes both IR-free theories and theories 
with weak IR fixed points. 

For QED3, the asymptotic freedom leads to fuv = f (4iV)-|-l [4]. The second 
term counts the single bosonic degree of freedom associated with a massless gauge 
boson in three space-time dimensions. In the first term, the AN counts the number 
of fermionic degrees of freedom associated with the N four-component fermions, 
and the 3/4 represents the Boltzman weighting of fermions in three space-time 
dimensions. 

The value of fm depends on whether the global symmetry is spontaneously 
broken. For large A'^, breaking is not expected, and only the massless fermions along 
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with the one gauge degree of freedom remain in the infrared spectrum. From Eq. 
2, one sees that the infrared theory is described by a weak (0(1/N)) fixed point. 
The free-field result, fiR = f(4Ar) + 1(= fuv), can then be corrected perturbatively 
in 1/N, the next-to-leading term being 0(1). The computation, similar in some 
ways to the corresponding perturbative computation in a 4D gauge theory, has not 
yet been done. The 4D result is negative, and if the same is true of the 3D 1/N 
computation, the inequality fiR < fuv will be satisfied. 

Now consider the more interesting possibility that for some finite N the global 
symmetry is broken. The fermions become massive and 2N'^ Goldstone bosons are 
formed. The infrared (massless) theory consists of only the Goldstone bosons and 
the gauge boson. Thus fjR = 2N^ + 1 and the inequality will be satisfied only 
if < 3/2. If this upper limit on Nc is correct, then QED3 cannot be in the 
broken phase if N is larger then |, meaning that the continuum gap equation [2] 
overestimates Nc- 

A natural question about this application of the inequality fjR < fuv has to 
do with the Mermin-Wagner-Coleman theorem [17]. We have determined fjR in the 
broken phase by computing the free energy at a temperature < T « a, counting 
the Nambu- Goldstone bosons plus the gauge field as the relevant, non-interacting 
degrees of freedom, multiplying by and an appropriate constant, and taking 

the limit T — > 0. Now at finite T, the IR behavior of this theory is that of the cor- 
responding 2D theory. But the Mermin-Wagner-Coleman theorem states that there 
can be no spontaneous symmetry breaking, with its Nambu-Goldstone bosons, in 
2D. This question was addressed by Rosenstein, Warr, and Park [18], who concluded 
that the symmetry is indeed unbroken at any non-zero T, and that for small T the 
zero-temperature Nambu-Goldstone bosons develop small masses. In our notation, 
one finds rrij^Q ~ T'^exp{—a/T), arising from the derivative interactions among the 
(pseudo) Nambu-Goldstone bosons. Since this mass vanishes more rapidly than T 
as r ^ 0, it has no effect on the computation of fiR. 

4 Lattice Studies 

Lattice studies may be the most reliable method to analyze dynamical mass genera- 
tion and directly determine Nc in QED3. Numerical simulations began over a decade 
ago, with a study in the quenched approximation giving preliminary evidence that 
the global chiral symmetry is spontaneously broken [19]. Recent advances in com- 
puting power have led to improved studies of mass generation in QED3. Simulations 
of 2-fiavor QED3 by Hands, Kogut and Strouthos [14] on lattices of up to 50^ sites 
report that the condensate is two orders of magnitude smaller than the quenched 
condensate value. They find that the value of the dimensionless condensate, defined 
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by 

^ = ^^, (4) 

where e is the dimensionful gauge couphng constant, is bounded above by 5 x 10~^. 
They also analyze finite-size effects and find that this bound is stable for lattice 
sizes ranging from 10'^ to 50^. They conclude that there is no decisive signal for 
chiral symmetry breaking for N >2. 

This conclusion is consistent with the conjectured inequality constraint de- 
scribed in the previous section, which led to N^. < 3/2 for QED3. On the other 
hand, this interpretation of the lattice results is likely premature. It has been noted 
by Gusynin and Reenders [20] that even a lattice of size 50^ provides a sufficient IR 
cutoff to affect substantially the value of Nc- They use the continuum gap equation 
first employed to note the existence of an Nc, modeling the effect of a finite lattice 
size by imposing an IR cutoff fi on the integral in the gap equation. Because of the 
scale-invariant form of the gap equation following from the dominance by the IR 
fixed point, the equation is logarithmic in character. A measure of the importance 
of the cutoff jU relative to the effective UV cutoff a is therefore In a/ fi. As this quan- 
tity decreases from infinity, a growing portion of momentum space is truncated, and 
therefore a stronger gauge coupling is required to break the symmetry. Thus, Nc 
should drop. Gusynin and Reenders find, for example, that with Ina/ ji still as large 
as 6 (a/jU 400), Nc drops by more than 30%. A lattice of size 50^ corresponds 
to a much smaller a///, indicating that chiral symmetry breaking should not take 
place for N = 2. 

If this analysis is qualitatively correct, that is, if the continuum IR cutoff 
correctly models the finite size of a lattice, then only a larger lattice would be able 
to determine accurately the value of Nc for a theory such as QED3, governed by 
a (conformal) IR fixed point. This same remark would apply to studies of the 
conformal phase transition as a function of the number of fermion species in four- 
dimensional gauge theories. 



5 Gap Equation Estimates 

Suppose that it becomes possible to carry out simulations with such large lattices 
that finite size effects are no longer important, and suppose that the results of the 
simulations for a (Eq. 4) with N = 2 continue to be bounded as in Ref. [14]. What 
would one conclude? Here we present a rough argument using the continuum gap 
equation indicating that even then it could be unclear whether iV = 2 is in the 
symmetric or broken phase. 

The spirit of the gap equation approach is to use the large-N form of the 
kernel of this equation (the large-N gauge-boson propagator ), and then to continue 



6 



to finite A^. The reliabifity of this approach, which first suggested the existence 
of a finite Nc 128/37r^), is not clear because the higher terms in the kernel are 
not parametrically small. There is, however, some evidence that the corrections are 
small numerically [3]. We use this approach to estimate the condensate for N = 2, 
expecting only that the estimate is order-of-magnitudc. 

The theory is UV-complete, rapidly damped at momentum scales beyond a. 
Thus, to a good approximation, the gap equation can be written down employing 
the a ^ oo form of the large-N kernel, with a then used as a UV cutoff. The 
resulting equation is 

where we make use of the non-local Kondo-Nakatani gauge for which the leading 
largc-N form of the wave function rcnormalization is unity. The critical N deter- 
mined from this equation is = 128/37r^. 

Numerical and analytical solutions of this equation [21] show that for a range 
of N such that ^ > 1, the overall scale of the solution of Eq.(5), set by S(0), 
is much less than a, and for A?^ = 2, 

^ 1.7 X 10-2. (6) 
a 

In the range of momentum S(p) < p < a the solution has the following approximate 
form 

EW = 5^fi»i„[i^^(M^) + *)l, (7) 

where S(0) is used to scale the log and 5 is a phase. 

We use these formulae to estimate the dimensionless fermion condensate of 
Eq. 4. The fermion condensate < V'V' > is given by the integral 

Noting that the sine function is bounded above by 1 in the region of integration up 
to a , we see that the dimensionless condensate a of Eq. 4 (=< N'^^ipil) /QAa^ >) is 
bounded above by 

167r2^ a ' ■ 

With N = 2,Nc = 128/37r2 and ^ ^ 1.7 x 10"^, this expression becomes 6 x 10"^ 
which is slightly above the upper bound of 5 x 10"^ found in lattice simulations [14]. 
We have numerically estimated the actual value of the condensate to be 4.4 x 10~^ 
which falls below the lattice bound [22] . 



7 



This suggests that a dimensionless condensate (4) of order 10~^ may nat- 
urally arise in QED3 for N of order 2 even when Nc is as large as 128/37r^ 4. 
It indicates that even if the numerical simulations are performed on a very large 
lattice, so that finite size effects are unimportant, a precise computation will be 
required to decide whether N = 2 is in the broken or symmetric phase. 



6 Abelian Higgs Model 

We conclude by discussing fermion mass generation when the Abelian gauge field 
has a non-zero mass. This theory has been used as a continuum description of the 
competition between long range antiferromagnetic order and superconducting order 
in planar cuprate systems [23], [24]. It can also be used to interpolate between 
QED3 and the 2-1-1 dimensional Thirring model, the latter having been studied in 
recent lattice simulations [14]. ( Related studies, using continuum gap equation 
methods and invoking the idea of hidden local gauge symmetry may be found in 
references [25] [26] ). Finally, it can be used to study the effect of an IR cutoff on 
QED3 [24] , similarly to the use of a simple IR cutoff as described in Section 3 [20] . 

The Lagrangian of the model is 

C = Cqedz + ^D^^^Di^^ - A($*$ - Nvf (9) 

where Cqed3 is given in equation (1), and $ is a complex scalar field with vacuum 
expectation value v. The addition of the scalar field does not change the global 
symmetry structure of the theory. We remove the Higgs boson from the spectrum 
by taking the limit A — > oo with v fixed. This leaves the theory UV-complete. The 
gauge boson mass is M = ey/Nv. To explore the phase structure as a function of 
N, we keep fixed the quantities a = and v'^. The gauge boson mass M is then 
also fixed. 

As in QED3, we define a dimensionless running gauge coupling a{k). The 
coefficient of g/^i, in the Euclidean gauge-boson propagator takes the form (k^ + 
+ ak)~^ to leading order in 1/N, where k is the magnitude of the Euclidean 
three-momentum and the last term arises from the N fermion loops. The k^ki, term 
of the propagator is dropped since it couples to a conserved current. Multiplying 
the propagator by the coupling and by k to make it dimensionless , we have 



e 



k 8 k 



"^^^^ ~ k^ + M^ + {e^N/8)k N (A;2 /a) + G'^ + k' ^^^^ 

where G = a/M"^. This effective coupling vanishes in both the UV and IR limits. 
It reaches a maximum no larger than 8/N, attainable only when M « a. As ^ 
increases, the maximum value decreases monotonically. Thus, for any M/a, the 
model is weakly coupled in the large-N limit at all momentum scales. 
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It is reasonable to expect that, as in QED3, a parity conserving mass is 
generated for N below some critical value Nc- Assuming an A^^ to exist, what 
can it depend on? Since the forces are strongly damped at scales beyond a, no 
additional cutoff is necessary to define the theory. Thus Nc can be a function of 
only the dimensionless ratio M/a. This dependence describes the boundary of a 
phase diagram for the theory. 

Consider first the limit M/a — > 0. On the one hand, this may be thought 
of as taking M — > with a fixed, giving QED3. As discussed earlier, estimates of 
Nc in this limit have ranged from 3/2 to approximately 4. The limit M/a 
may also be taken by letting a — > oo with G = a/M^ fixed. This leads to the 2+1 
dimensional Thirring model, meaning that Nc for this model is the same as that for 
QED3. This assumes that the Thirring model is treated with any UV cutoff taken 
very large relative to G~^. For finite N, this can be implemented using lattice 
techniques. (For large N, these two ways of taking the limit M/a — >■ may be 
considered by examining a{k) (Eq. 11). In the case M — > with a fixed (QED3), 
there appears an IR {k « a) fixed point of strength 8/N. In the case a ^ oo with 
fixed , there is a UV (A; >> G"^) fixed point of the same strength). 

As M/a is increased from zero, low momentum components are damped out. 
If the force driving the symmetry breaking is attractive at all scales (as in the large 
N approximation), a stronger infrared coupling is required to trigger dynamical 
symmetry breaking, and therefore Nc will decrease. At some value of M/a, Nc 
drops below unity meaning that symmetry breaking will not take place at all [23]. 
This critical curve should be similar to that determined by Gusynin and Reenders 
[20], with the gauge boson mass replacing the explicit IR cutoff. The character of 
the phase transition along this critical curve is also of interest. For M/a = 0, as 
discussed earlier, because the force is of infinite range, the continuum gap equation 
suggests that the transition is of infinite order [13]. With M > 0, the force becomes 
of finite range, and the transition is of second order. 

Does the conjectured inequality fjR < fjjv provide any information about Nc 
as a function of M/al With a = finite, the theory is asymptotically free. A free- 
field computation then gives fuv = |(4A") + 2 = 3A^ + 2 [4]. The second term counts 
the 2 bosonic degrees of freedom associated with a massive gauge boson in three 
space-time dimensions. As in QED3, the value of fm depends on whether the global 
symmetry is spontaneously broken. For large N, breaking is not expected, and, if 
M is nonzero, the infrared theory consists of only the (non-interacting) massless 
fermions. Then fjji = |(4A^), and the inequality is satisfied. Now suppose that for 
some finite N the global symmetry is broken. The fermions become massive, and for 
finite gauge boson mass M the IR theory consists of only the 2A?"^ non-interacting 
Goldstone degrees of freedom. Thus fiR = 2A^^, and the inequality demands that 
Nc be no greater than 2. But the inequality already demands that Nc be no greater 
than 3/2 when M = 0, and we expect that Nc will decrease as M/a increases. So 
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it would seem that the conjectured inequahty has nothing useful to say about the 
shape of this critical curve. 

Lattice simulations of the Abelian Higgs model for an appropriate range of 
values of M/a could determine the shape of this critical curve. Simulations have 
been done only for QED3 and the Thirring model, which correspond to the limit 
M/a — > 0. We discussed in Section 3 the sensitivity of the QED3 simulation to finite 
size effects. A numerical study of the Thrring model was reported by Hands and 
Lucini [27]. We noted above that providing this model is treated with a UV cutoff 
large compared to the inverse four-fermion coupling , the value of the critical 
coupling Nc should approach that of QED3. The lattice simulations of Ref. [27] do 
take the ultraviolet cutoff larger than G~^, and find chiral symmetry breaking for 
iV Ri 3. Their results remain sensitive to the UV cutoff, but we anticipate that Nc 
will only increase as the cutoff is made larger (as the continuum limit is approached). 
This suggests that Nc > 3 for the continuum Thirring model, and therefore also for 
QED3. 

We end this section by digressing to use the large-N Thirring model to explain 
why the inequality [4] was restricted to asymptotically free theories. The model 
becomes non-interacting in the infrared (A; << G~^) so fjji = 3N. Having taken 
a — > GO, the resultant weak UV fixed point means that fuv exists and may be 
computed in the expansion. To leading order, fuv = 3A?^. The computation of 
the next-to-leading corrections should be the same as that of fm in QED3, since the 
latter is governed by an IR fixed point identical in strength to the UV fixed point 
of the Thirring model. If the result is negative as anticipated, then fjR > fuv 
for this model. Another 2-1-1 dimensional model with a UV fixed point, leading to 
fiR > fuv, was studied by Sachdev [28]. In light of these examples, the conjectured 
inequalty was restricted to asymptotically free theories. 

7 Summary and Discussion 

A key question for a non-compact Abelian gauge theory in three space-time dimen- 
sions is the critical number Nc of four-component fermions below which there is 
fermion mass generation with an associated chiral symmetry breaking, and above 
which the fermions remain massless. Over the years, estimates of Nc with the 
gauge symmetry unbroken (QED3) have ranged from infinity to one. Recent stud- 
ies agree that it is finite and about 4, but some estimates continue to be as low as 
1. Since this theory is UV-complete, damping rapidly at momentum scales above 
a = e^N/S, the value of Nc depends only of the IR behavior of the theory, where 
conformal symmetry sets in. 

This question has a potential real-world interest, because of the possibility 
that QED3 with N = 2 can describe the physics of the planar anti-ferromagnet [5] 
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[6]. The observed antiferromagnetic (Neel) ordering eorresponds to chiral symmetry 
breaking in QED3 [7] [8] . Numerical simulations [29] of the planar antiferromagnet 
also indicate Neel ordering. Herbut and collaborators [12] have challenged the idea 
that non-compact QED3 is relevant to the planar antiferromagnet, noting that the 
underlying theory is compact QED3 with its topological (instanton) structure. Some 
argue, however, that for large enough N, instantons and anti-instantons are bound 
at distance scales of order \/a by a logarithmic term in the action [10]. Then for 
N = 2 the non-compact theory could be employed at larger distances. But Herbut 
et al [12] claim that the finite density of instantons and anti-instantons screen the 
logarithmic term in the action and that instantons axe therefore unsuppressed at 
all N. If this is the case, the non-compact theory cannot be used for this condensed 
matter system [30]. 

In general, the question of whether the value of Nc as determined by (UV- 
complete) QED3 is relevant to condensed matter systems depends on the role of 
new, short-distance physics. With the physical lattice spacing of order 1/a = N/Se^, 
physics at this scale can be represented in the effective theory at momentum scales 
below a by the addition of a set of higher-dimension operators. The importance 
of these terms depends on their strength. We have analyzed this question in the 
framework of the continuum gap equation in 4-dimensions [31], which has much the 
same structure as QED3 treated in the 1/N expansion. With a four-fermion additive 
term, there is a critical four-fermion coupling (attractive) of order unity above which 
this term by itself can trigger spontaneous chiral symmetry breaking. As long as the 
magnitude of this coupling, either attractive or repulsive, is less than one quarter 
of its critical value [32, 33], the additive term is unimportant, not affecting the 
value of Nc- This limit may not be quantitatively accurate, but it seems reasonable 
that a limit of this order exists. Whether this limit is satisfied for condensed matter 
systems such as the planar antiferromagnet is not clear. Thus even if the non- 
compact theory does provide a continuum descriptiion of the antiferromagnet, the 
value of Nc for this system may or may not be determined accurately by QED3 
itself. 

Returning to our summary of QED3, we have briefly reviewed the various 
approaches to the determination of N^. A recently conjectured inequality constrain- 
ing asymptotically free quantum field theories [4], when applied to QED3, suggests 
that Nc < 3/2. This bound disagrees with the estimate A^c ~ 4 emerging from 
the continuum gap equation approach which originally suggested the existence of a 
finite Nc- 

Recent lattice simulations of A^ = 2 QED3 [14] indicate that the fermion 
condensate is very small, possibly signaling that QED3 is in the symmetric phase, 
and therefore that Nc < 2. This conclusion depends, however, on the importance 
of finite-size effects. A recent analysis of Gusynin and Reenders [20] using the 
continuum gap equation indicates that a lattice size of 50^, used in the simulations. 
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provides a strong enough suppression of IR physics to reduce Nc below 2, thus 
accounting for the obseravation of a small upper bound for the condensate. Using 
the solution of the gap equation with no IR cutoff, we have made a rough estimate of 
the fcrmion condensate, assuming that Nc is larger - approximately 4 as suggested 
by the continuum gap equation studies [3] . Interestingly, even though A'' = 2 is then 
well into the broken phase, the condensate is estimated to be very small, below the 
upper bound of Ref. [14] . This indicates that even if simulations are conducted on a 
large enough lattice to make finite size effects unimportant, a precise determination 
of the condensate is needed to decide whether QED3 with iV = 2 is in the symmetric 
or broken phase. 

The value of Nc in the Abelian Higgs model is also of interest. This model 
leads to either QED3 or the 2+1-dimensional Thirring model in the limit M/a — ^ 
depending on how the limit is taken. But since Nc depends on only the ratio M/a 
the value of Nc must be the same in both models. Lattice simulations of the Thirring 
model on a 12^ lattice lead to Nc^ 3 [27]. This value remains sensitive to the UV 
cutoff (the lattice spacing), and should increase toward the Nc of QED3 as the 
continuum limit is approached. On the other hand the inequality (3) applied to the 
Abelian Higgs model yields Nc <2. 

To summarize, after more than fifteen years there is still no definitive answer 
to the question of the critical number of four-component fermions in QED3 (and its 
generalization to the Abelian Higgs model) marking the boundary between broken 
and unbroken chiral symmetry. This question is of interest because it could be 
relevant to the behavior of certain condensed matter systems and because the answer 
requires an understanding of a strongly coupled, UV-complete quantum field theory. 

8 Acknowledgements 

The work described here was supported in part by funds provided by the U.S. De- 
partment of Energy under Contracts DE-AC02-76ER0 3075 and DE-FG02-84ER40153. 
We thank Simon Hands and Subir Sachdev for valuable conversations, Frank Pinski 
for his assistance with computing and V. P. Gusynin for helpful correspondence. 

References 

[1] R.Pisarski, Phys. Rev. D29. 2423 (1984). 

[2] Thomas Appelquist, Daniel Nash and L.C.R. Wijewardhana , Phys. Rev. Lett. 
60, 2575 (1988). 

[3] D. Nash, Phys. Lett. 62, 3024 (1989). 



12 



[4] T. Appclquist, A. Cohen, and M. Schmaltz, Phys. Rev. D60, 045003 (1999). 

[5] I. Affleck and J. B. Marston, Phys. Rev B39, 11538 (1989), L.B.Ioffe and A. I. 
Larkin, Phys. Rev. B39, 8988 (1989). 

[6] M.C.Diamantini, P.Sodano and G.Semenoff, Nuc.Phys.B406,595(1993) . 

[7] J. B. Marston, Phys. Rev. Lett. 64, 1166 (1990). 

[8] R.Laughhn and Z.Zou, Phys. Rev. B41, 664 (1990). 

[9] D. H. Kim and Patrick A. Lee, cond-matt/9810130 pubhshed in Ann. Phys. 
(N. Y.) 272, 130 (1999). 

[10] H. Kleinert, F. Nogueira, and A. Sudbo, Phys. Rev. Lett. 88 (2002) 232001, 
hep-th/0201168. 

[11] Z. Tesanovic, O. Vafek and M. Franz, Phys. Rev. Lett. 87, 257003 (2001) , 
Phys. Rev. B 65, 180511 (2002), M. Franz, Z. Tesanovic and O. Vafek, Phys. 
Rev. B 66, 054535 (2002), L F. Herbut, Phys. Rev. B66, 094504 (2002). 

[12] I. Herbut and B. Seradjeh, Phys. Rev. Lett. 91, 171601, 2003 , cond- 
mat/0305296, M. Case, B. Seradjeh, and L Herbut, Nucl. Phys. B676: 572-586, 
2004, e-print archive cond-mat/0308260. 

[13] T. Appelquist, J. Terning, and L. C. R. Wijewardhana, Phys. Rev. Lett. 75, 
2081 (1995). 

[14] S. Hands, J. Kogut and K. Strouthos, Nucl. Phys. B645: 321-336, 2002; e-Print 
Archive: hep-lat/0208030. 

[15] D. Lee, P. Maris, Phys. Rev. D67: 076002, 2003, e-print archive hep- 
lat/0212033. 

[16] C. Vafa and E. Witten, Commun. Math. Phys. 95, 257 (1984); Phys. Rev. 
Lett. 53, 535 (1984) and T. Appelquist, M. Bowick, D. Karabali and L. C. R. 
Wijewardhana, Phys. Rev. D33, 3704, 1986. 

[17] By N.D. Mermin, H. Wagner, Phys. Rev. Lett. 17: 1133-1136,1966, S.Coleman, 
Commun. Math. Phys.31 :259-264, 1973. 

[18] B. Rosenstein, B. Warr, and S. Park, Nuclear Physics B336, 435 (1990). 

[19] S. Hands and J. Kogut, Nuc. Phys. B335 (1990) 455. 

[20] V.P.Gusynin and M.Reenders, Phys. Rev. D68, 025017 (2003) , hep- 
ph/0304302. 



13 



[21] Thomas Appelquist, Mark J. Bowick, Eugene Cohler, L.C.R. Wijewardhana 
Phys. Rev. Lett. 55:1715, 1985, Thomas W. Appelquist, Mark J. Bowick, Dim- 
itra Karabah, L.C.R. Wijewardhana Phys. Rev. D33, 3704,1986. 

[22] V.Gusynin has also obtained a similar estimate for the two flavor condensate, 
using analytical solutions to the gap equation described in Ref. [20]. 

[23] Guo-Zhu Liu and Ceng Cheng, Phys. Rev. D67, 065010, 2003 e-Print Archive: 
hep-th/0211231. 

[24] T. Pereg-Barnea and M. Franz, Phys. Rev. B67 060503 (2003), cond- 
mat/0209301. 

[25] Taichi Itoh, Yoonbai Kim, Masaki Sugiura, Koichi Yamawaki, Prog. Theor. 
Phys. 93: 417-440, 1995, hep-th/9411201. 

[26] M.Sugiura, Prog. Theor. Phys. 97, 311(1997), hcp-th/9611198. 

[27] Simon Hands, Biagio Lucini, Phys. Lett. B461 :263-269,1999, e-Print Archive: 
hep-lat/9906008. 

[28] Subir Sachdev, Phys. Lett. B309, 285-288, 1993, e-Print Archive: hep- 
th/9305131. 

[29] U.-J. Wiese and H.-P. Ying, cond-mat/9212006 and B. B. Beard and U.-J. 
Wiese, Phys. Rev. Lett. 77 , 5130 (1996). 

[30] Some discussion of these issues may be found in T. Senthil, L. Balents, S. 
Sachdev, A. Vishwanath, and M. P. A. Fisher, cond-mat/0312617. 

[31] T. Appelquist, J. Terning, and L. C. R. Wijewardhana, Phys. Rev. D44, 871 
(1991). 

[32] T.Appelauist, M.Soldate, T.Takeuchi and L.C.R.Wijewardhana, Proceedings 
of the Johns Hopkins workshop on current problems in physics, edited by 
G.Domokos and S.Kovesi-Domokos, World Scientific, Singapore 1988, Kei-ichi 
Kondo , Hidetoshi Mino, Koichi Yamawaki, Phys. Rev. D39, 2430,1989. 

[33] M. Carena, T.E.Clark and C.E.M. Wagner Nucl. Phys. B356,117-153,1991, U. 
Mahanta Phys. Rev. D44, 3356-3359, 1991. 



14 



